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Abstract. We study systematic errors of the ground-state parameters obtained by Shifman- 
Vainshtein-Zakharov (SVZ) sum rules, making use of the harmonic-oscillator potential model as 
an example. In this case, one knows the exact solution for the polarization operator, which allows 
one to obtain both the OPE to any order and the parameters (masses and decay constants) of the 
bound states. We determine the parameters of the ground state making use of the standard proce- 
dures of the method of sum rules, and compare the obtained results with the known exact values. We 
' show that in the situation when the continuum contribution to the polarization operator is not known 

and is modelled by an effective continuum, the method of sum rules does not allow to control the 
systematic uncertainties of the extracted ground-state parameters. 
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> ' A QCD sum-rule calculation of hadron parameters [1] involves two steps: (i) one 
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calculates the operator product expansion (OPE) series for a relevant correlator, and (ii) 
one extracts the parameters of the ground state by a numerical procedure. Each of these 
steps leads to uncertainties in the final result. 

■ The first step lies fully within QCD and allows a rigorous treatment of the uncertain- 
O I ties: the correlator in QCD is not known precisely (because of uncertainties in quark 
^ I masses, condensates, a^, radiative corrections, etc), but the corresponding errors in the 

correlator may be systematically controlled (at least in principle). 

The second step lies beyond QCD and is more cumbersome: even if several terms 
^ I of the OPE for the correlator were known precisely, the hadronic parameters might be 

■ extracted by a sum rule only within some error, which may be treated as a systematic 
error of the method. It is useful to recall that a successful extraction of the hadronic 
parameters by a sum rule is not guaranteed: as noticed already in the classical papers 
[1,2], the method may work in some cases and fail in others; moreover, error estimates 
(in the mathematical sense) for the numbers obtained by sum rules may not be easily 
provided — e.g., according to [1], any value obtained by varying the parameters in the 
sum-rule stability region has equal probability. However, for many applications of sum 
rules, especially in flavor physics, one needs rigorous error estimates of the theoretical 
results for comparing theoretical predictions with the experimental data. Systematic 
errors of the sum-rule results are usually estimated by varying the Borel parameter and 
the continuum threshold within some ranges and are believed to be under control. 

Here we present the results of our recent analysis of the systematic uncertainties of the 
sum-rule procedure [3]. To this end, a quantum- mechanical harmonic-oscillator (HO) 



potential model is a perfect tool: in this model both the spectrum of bound states (masses 
and wave functions) and the exact correlator (and hence its OPE to any order) are known 
precisely. Therefore, one may apply the sum-rule machinery for extracting parameters 
of the ground state and test the accuracy of the extracted values by comparing with the 
exact known results. In this way the accuracy of the method can be probed. 

The model. To illustrate the essential features of the QCD calculation, let us consider a 
non-relativistic model with a potential containing both a Coulomb and a confining part: 

H = Ho + V{r), HQ = f/2m, V(r) = a./r + W (1) 

The polarization operator n(£') is defined by the full Green function G{E) — {H — E)~^ 
as follows: 

Il{E) = {iK/mf^ {ff = Q\G{E)\ri = 0). (2) 

The full Green function satisfies the Lippmann-Schwinger operator equation 

G-\E)-G^\E)^V, (3) 

which may be solved by constructing the expansion in powers of the interaction V: 

G{E) = Go{E)-Go{E)VGo{E) + ---. (4) 

Respectively, for the polarization operator n(£') one may construct the expansion shown 
in Fig. 1. If the confining potential is known, one can calculate each term of the expan- 
sion. If one does not know the confining potential precisely, an explicit calculation is not 
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possible. In this case one can explicitly calculate the contribution IIq + 11" + 11" H 

of those diagrams which do not contain the confining potential (first line in Fig. 1) and 
parametrize other diagrams, in which the confining potential appears, taking into ac- 
count the symmetries of the theory. In QCD such diagrams are parametrized in terms of 
condensates and radiative corrections to them. 

In what follows we make a further simplification of the model and switch off the 
Coulomb potential, which we, in principle, know how to deal with. We shall retain only 
the confining potential, choose it as the HO potential 

V{r) = m(0^r^/2, r= |r|, (5) 

and study n(ju) [the Borel transform of the polarization operator !!(£')] which gives the 
evolution operator in the imaginary time 1 /jU : 

n(;u) = (2;r/m)3/2 ^ q| exp(-///;u)|r,- = 0). (6) 



For the HO potential (5), the exact analytic expression for n(ju) is known [2]: 
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FIGURE 1. The expansion of the polarization operator n(/?^) in powers of the interaction: wavy Hnes 
(green) represent the Coulomb-part exchange, dashed lines (red) represent the confining-part exchange. 



Expanding this expression in inverse powers of [l, we get the OPE series for Tl{[i): 



nopE(ju) = no(ju) +ni (ju) + n2(ju) + ■ ■ ■ = ju 
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and higher coefficients may be obtained from (7). Each term of this expansion may be 
also calculated from (2) and (4), with Do corresponding to Gq. 

The "phenomenological" representation for n(jU) is obtained by using the basis of 
hadron eigenstates of the model, namely 

oo 

n(Ai) = £i?„exp(-£„/Ai), (9) 

)j=0 

with E„ the energy of the n-th bound state and R„ (the square of the leptonic decay 
constant of the n-th bound state) given by 

i?„ = (2;r/m)3/2|vi/„(F = o)|2. (10) 

For the lowest states, one finds from (7) 

3 7 
Eq = -(0,Rq = 2^2(0^/^, £i = -(0, i?i = 3720)3/2, .... (H) 

Now, imagine we know n(/i) numerically and want to extract the parameters of the 
ground state Eq and Rq. At large values of the Euclidean time l/ji the polarization op- 
erator is dominated by the ground state. Therefore, one has a plateau in — ^logn(ju) 
and in the product exp(£'o//i)n(/i) (Fig. 2). In practice, working at sufficiently large 
values of l//i, one constructs — ^ logn(/i) and, after verifying that the plateau has al- 
ready been reached, determines Eq. Then one studies n(/i) exp(£'o/M) and reads off Rq. 
In principle, it is possible to reach the plateau if one retains sufficiently many terms of 
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FIGURE 2. Dashed (black) line - the calculation for the exact n(/i); dotted (red) line - the calculation 
making use of the OPE for n(/x) with 4 power corrections, solid (green) line - the same OPE with 100 
power corrections. 



the truncated OPE: the radius of convergence of the truncated OPE to the exact function 
n(/i) increases with the increase of the number of terms taken into account, and, e.g., 
for > 100 the range of convergence already reaches the plateau (Fig. 2). 

In practice, however, one can calculate only the first few terms and, consequently, the 
radius of convergence is rather small. For instance, with 4 terms the convergence region 
extends only up to values of l//i for which the ground state gives 85% of the correlator. 

The method of SVZ sum rules aims at extracting the parameters of the ground state 
from the first few terms of the OPE, e.g. in the region where the contribution of higher 
states is by far not negligible. 

Sum rule. The sum rule claims the equality of the correlator calculated in the "quark" 
basis and in the hadron basis: 
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rfzpo(z)exp(-z//i) +11^^^ 
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(12) 



Following [2], we use explicit expressions for the power corrections, but for the zero- 
order free-particle term we use its expression in terms of the spectral integral. The reason 
for this will become clear in a few lines. 

Let us introduce the effective continuum threshold Zeff(M)j different from the physical 
jLi -independent continuum threshold Zcont, by the relation 

oo oo 

ncont(M) = J dzppheniz)exp{-z/{i)= J Jzpo(z) exp(-z/^). (13) 



^cont 



Zeff(M) 



Generally speaking, the spectral densities Pphen(z) and po{z) are different functions, so 
the two sides of (13) may be equal to each other only if the effective continuum threshold 



depends on /i. In our model, we can calculate Hcont precisely, as the difference between 
the known exact correlator and the known ground-state contribution, and therefore 
we can obtain the function Zeff(M) by solving (13). In the general case of a sum-rule 
analysis, the effective continuum threshold is not known and is one of the essential 
fitting parameters. 

Making use of (13), we rewrite now the sum rule (12) in the form 

i?oexp(-£'o/M) = n(jU,Zeff(M)), (14) 
where the cut correlator n(/i,Zeff(ju)) reads 

Zeff(M) 

h-^/verxnf — 7 / II.) -\- 11^/^ - ^ , 

4;u2 480 

u 

As is obvious from (14), the cut correlator satisfies the equation 
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n(jU,Zeff(M)) = -7^ J dzVzew{-z/ H) + H^^^ 



(15) 



logn(Ai,Zeff(ju))=^o. (16) 



The cut correlator is the actual quantity which governs the extraction of the ground-state 
parameters. 

The sum rule (14) allows us to restrict the structure of the effective continuum 
threshold Zeff(M)- Let us expand both sides of (14) near co/n = 0. The l.h.s. contains 
only even powers of y^co/Ji; power corrections on the r.h.s. contain only odd powers 
of ^/(o/J^. In order that both sides match each other, the effective continuum threshold 
cannot be constant but should be a power series of the parameter y^co/fi: 
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(17) 



Inserting this series in (14) and expanding the integral on the r.h.s., we obtain an infinite 
chain of equations emerging at different orders of ^J~&J\i. These equations allow us 
to obtain for any £0 and i?o within a broad range of values a solution Zeif{^-,RoiEo) 
which exactly solves the sum rule (12). Therefore, in a limited range of ju the OPE 
alone cannot say much about the ground-state parameters. What really matters is the 
continuum contribution, or, equivalently, Zeff(M)- Without constraints on the effective 
continuum threshold the results obtained from the OPE are not restrictive. 

The approximate extraction of Eq and Rq worked out in a limited range of values 
of jU becomes possible only by constraining Zesin). If the constraints are realistic and 
turn out to reproduce with a reasonable accuracy the exact ZefflM)^ then the approximate 
procedure works well. If a good approximation is not found, the approximate procedure 
fails to reproduce the true value. Anyway, the accuracy of the extracted value is difficult 
to be kept under control. This conclusion is quite different from the results of QCD sum 
rules presented in the literature (see e.g. the review [4]). We shall demonstrate that a 
typical sum-rule analysis contains additional explicit or implicit assumptions and criteria 
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FIGURE 3. Different choices of the effective continuum threshold Zeff(M) (a) and the corresponding 
R{fi) (b): 1 [soHd (red) Hne] the exact effective continuum threshold as obtained by a numerical solution 
of (18) for E ^ Eq and R = Rq. 2 [long-dashed (blue) line] the effective continuum threshold obtained by 
solving the sum rule (18) for E ~ Eq and R — QJRq, 3 [dash-dotted (green) line] same as line 2, but for 
for£' = £'oand/?=1.15/?o. 



for extracting the parameters of the ground state. Whereas these assumptions may lead 
to good central values of hadron parameters, the accuracy of the extracted values cannot 
be controlled within the method of QCD sum rules. 

Numerical analysis. In practice, one knows only the first few terms of the OPE, so one 
must stay in a region of /i bounded from below to guarantee that the truncated OPE series 
reproduces the exact correlator within a controlled accuracy. The "fiducial" [1] range of 
jl is the range where, on the one hand, the OPE reproduces the exact expression better 
than some given accuracy (e.g., within 0.5%) and, on the other hand, the ground state 
is expected to give a sizable contribution to the correlator. If we include the first three 
power corrections, 111, 112, and Ut,, then the fiducial region lies at (y//i < 1.2 (see Fig. 3). 
Since we know the ground-state parameters, we fix (o//i > 0.7, where the ground state 
gives more than 60% of the full correlator. So our working range is 0.7 < (D//i < 1.2. 

Obviously, if one knows the continuum contribution with a reasonable accuracy, one 
can extract the resonance parameters from the sum rule (12). We shall be interested, 
however, in the situation when the hadron continuum is not known, which is a typical 
situation in heavy-hadron physics and in studying properties of exotic hadrons. Can we 
still extract the ground-state parameters? 

We shall seek the (approximate) solution to the equation 

oo 

Rexp{-E/{i)+ J Jzpo(z)exp(-z/Ai) =noPE(M) (18) 

Zeff(M) 

in the range 0.7 < co/jl < 1.2. Hereafter, we denote by E and R the values of the 
ground-state parameters as extracted from the sum rule (18). The notations Eq and Rq 
are reserved for the known exact values. 

^-dependent ejfective continuum threshold. In many interesting cases the ground-state 
energy can be determined, for instance, from experiment. However, fixing the ground- 
state energy E equal to its known value Eq does not help: for any R within a broad range 
one can still find a solution Zeff{lJ-,R) which solves the sum rule (18) exactly (Fig. 3). In 
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FIGURE 4. Constant effective continuum threshold Zc- for three different values of Zc (a) and the 
corresponding (b). 



order to extract R, one needs to constrain the effective continuum threshold; moreover, 
this constraint determines the actual value of R which one obtains from the sum rule. If 
one requires, e.g., Zeff(M) > for all /i, then the sum rule (18) may be solved for any 
R within the broad range 0.7 < R/Rq < 1. 15 (Fig. 3). Thus, the sum rule alone, without 
knowing the continuum contribution, cannot determine R. 

Constant effective continuum threshold. Strictly speaking, a constant effective contin- 
uum threshold Zeff (m) = Zc = const is incompatible with the sum rule. Nevertheless, this 
Ansatz may work well, especially in our model: as can be seen from Fig. 3(a), the exact 
Zeff (m) is almost flat in the fiducial interval. Therefore, the HO model represents a very 
favorable situation for applying the QCD sum-rule machinery. 

Now, one needs to impose a criterion for fixing Zc- A widely used procedure is the 
following [5]: one calculates 
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(19) 



which now depends on /i due to approximating Zeff(M) by some constant. Then, one 
determines /io and Zc as the solution to the system of equations 
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(20) 



yielding Zc = 2.454 ft), /io/ft) = 1 (Fig. 4), and a good central-value estimate R/Rq = 

0. 96, with R being extremely stable in the fiducial range. 

Note, however, a dangerous point: (i) a perfect description of n(ju) with better than 
1% accuracy, (ii) the deviation of the E{ii,Zc) from Eq at the level of only 1%, (iii) and 
extreme stability of i?(/i) in the full fiducial range leads to a 4% error in the extracted 
value of R\ Clearly, this error could not be guessed on the basis of the other numbers 
obtained, and it would be wrong to estimate the error, e.g., from the range covered by R 
when varying the Borel parameter /i within the fiducial interval. 

Let us summarize the lessons we have learnt: 

1. The knowledge of the correlator in a limited range of the Borel parameter ji is not 
sufficient for an extraction of the ground- state parameters with a controlled accuracy: 



Rather different models for the correlator in the form of a ground state plus an effective 
continuum lead to the same correlator. 

2. The procedure of fixing the effective continuum threshold by requiring the average 
mass calculated with the cut correlator to reproduce the known ground-state mass [5-7] 
is, not restrictive: a /i-dependent effective continuum threshold Zeff(M) which solves the 
sum rule (18) leads to the cut correlator (15) which automatically (i) reproduces exactly 

= £0 for all values of the Borel parameter ji, and (ii) yields a -independent value 
of R which, however, may be rather far from the true value. 

3. Modeling the hadron continuum by a constant effective continuum threshold allows 
one to determine Zc by, e.g., requiring the average energy E(ii) to be close to Eq in 
the stability region. In the model under discussion (where the exact effective continuum 
threshold is in fact almost constant) one obtains in this way a good estimate i?/i?o = 0.96, 
with practically /i -independent R. The unpleasant feature is that the deviation of R from 
Rq turns out to be much larger than the variations of E{ii) and R over the fiducial range. 
Therefore, the systematic errors cannot be estimated by looking at the variation of R 
within the fiducial domain of ji. 

Thus, we conclude that in those cases where the continuum contribution is not known 
(which is a typical situation encountered, e.g., in analyses of heavy-meson observables) 
the standard procedures adopted in a QCD sum-rule extraction of hadron parameters do 
not allow to control the systematic uncertainties. Consequently, no systematic errors for 
hadron parameters obtained with sum rules can be provided. Let us emphasize that the 
independence of the extracted values of hadron parameters from the Borel mass does not 
guarantee the extraction of their true values. 

Nevertheless, in the model under consideration the sum rules give good estimates for 
the parameter Rq. This may be due to the following specific features of the model: (i) 
a large gap between the ground state and the first excitation that contributes to the sum 
rule; (ii) an almost constant exact effective continuum threshold in a wide range of jU. 
Whether or not the same good accuracy may be achieved in QCD, where the features 
mentioned above are absent, is not obvious, and requires more detailed investigations. 

This shortcoming — the impossibility to control the systematic errors — remains the 
weak feature of the method of sum rules and an obstacle for using the results from QCD 
sum rules for precision physics, such as electroweak physics. 
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